در این تحقیق، پدیده خزش دیسک دوار پلیمری با استفاده از مدل ویسکو-هایپرالاستیک ماکسول تعمیم یافته مورد بررسی قرار گرفته است. پس از استخراج معادله دیفرانسیل جزئی تعادل لاگرانژی حاکم بر مسئله، تحلیل دیسک دوار با کدنویسی در FlexPDE انجام شد. مدل سازی دیسک در ANSYS با کدنویسی در محیط APDL، نشان داد که جابجایی شعاعی و تنش ون-میزز، تطابق بسیار خوبی با نتایج FlexPDE دارند. مزایای کدنویسی FlexPDE نسبت به انسیس از جمله تحلیل یک بعدی تنش-صفحه ای محورمتقارن بجای تحلیل دوبعدی، کاهش هزینه محاسباتی، امکان تعریف ضخامت متغیر (بدون کدنویسی اضافی) و نیاز به تعداد المان های کمتر در راستای شعاعی برای رسیدن به دقت مناسب (لزوم بکارگیری 20 المان در FlexPDE در مقایسه با 100 المان در ANSYS)، می باشد.نشان داده شد با گذر زمان و افزایش سرعت زاویه ای، جابجایی شعاعی و تنش ون-میزز دیسک دوار ناشی از پدیدۀ خزش، افزایش می یابد. مشخص شد که با افزایش سرعت زاویه ای و کاهش توان در پروفیل ضخامت n_h، مقدار جابجایی و تنش ون-میزز در یک زمان مشخص، افزایش می یابد؛ اما تغییر در سرعت زاویه ای (به عنوان بار اعمالی) و تغییر در پارامتر n_h (به عنوان یک مشخصه ی هندسی)، تأثیر چندانی در زمان آرامش دیسک دوار ندارد.